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The transition to turbulence in a precessing cylindrical vessel is experimentally inves¬ 
tigated. Our measurements are performed for a nearly-resonant configuration with 
an initially laminar flow dominated by an inertial mode with azimuthal wave number 
m = 1 superimposed on a solid body rotation. By increasing the precession ratio, we 
observe a transition from the laminar to a non-linear regime, which then breakdowns 
to turbulence for larger precession ratio. Our measurements show that the transi¬ 
tion to turbulence is subcritical, with a discontinuity of the wall-pressure and the 
power consumption at the threshold eLT- The turbulence is self-sustained below this 
threshold, describing a bifurcation diagram with a hysteresis. In this range of the 
control parameters, the turbulent flows can suddenly collapse after a hnite duration, 
leading to a dehnitive relaminarization of the flow. The average lifetime (r) of the 
turbulence increases rapidly when e tends to elt- 
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I. INTRODUCTION 


Rotating flows are ubiquitous in geophysical, astrophysical and industrial context. When 
an external torque is applied perpendicularly to the rotation axis of a solid body, it displays 
a processing motion in order to conserve the angular momentum. The precession refers to 
the continuous change of direction of the axis of rotation, which rotates around the axis of 
precession. The response of a fluid flow driven by the precession of the vessel is a complex 
process^, which results mostly in the interplay of inertial waves, Ekman boundary layers 
and base flow. It has been observed experimentally that the fluid flow driven by precession 
can amplify^ a magnetic held, and simulations have demonstrated that a dynamo effect can 
operate^i^. The present study belongs to the framework of the project DRESDYN-^-, which 
aims to investigate experimentally a precession driven dynamo. 

The observed regimes and instabilities of a processing how depend strongly on the con- 
hguration, which includes the geometry (spherical, ellipsoidal or cylindrical), the rotation of 
the body, and the strength of the precession. For sake of clarity, we recall only the results 
in cylindrical geometry. In most of the experimental investigation, the set-up is composed 
of a cylinder rotating around its revolution axis at the rotation rate flc and this cylinder is 
mounted on a turntable rotating at the precession rate Qf When the angle between the axes 
of rotation and precession is non-zero, the how deviates from the solid body rotation, due 
to the acceleration of the rotation axis. Initially observed by McEwaii^, the laminar how is 
composed of Kelvin modes superimposed on a solid body rotatioii^i^— . The Kelvin modes^ 
(or inertial modes) are the eigen-modes of a rotating cylindrical vessel. Each mode has its 
own eigen-frequency u, which depends on the radial, axial and azimuthal wave numbers. 
The gyroscopic force due to the precession excites the modes with azimuthal wave number 
m = 1 and angular frequency fic- A mode can be resonant if its eigen-frequency is close to 
the angular forcing frequency flc- Out of the resonance, the amplitude of the modes scales 
like \u: — uud it increases with the strength of the precession, quantihed by the pre¬ 

cession ratio e = VLt/Vtc- At the resonance, viscous and non-linear ehects play an important 
role in the saturation of the amplitude of the Kelvin modesl^'^i'^^'l^. 

Suggested by McEwan^, Kerswell^ has demonstrated that the Kelvin modes forced by 
precession become unstable by a parametric instability. The mechanism is based on a 
triadic resonance with two other free Kelvin modes, i.e. two modes not directly forced by 
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the precession. Lagrange et al.—^ confirmed experimentally this scenario for a processing 
flow with a single forced Kelvin mode close to its resonance. The parametric instability can 
either saturate^^ or trigger a secondary instability called resonant collapse^'i^.'^S. A weakly 
non-linear analysis performed by Lagrange et alM- predicts a subcritical bifurcation leading 
to a stationary state or a cyclic dynamics for larger Reynolds number. When the instability 
does not saturate, a sequence of bifurcations called resonant collapse triggers a breakdown 
of the flow into a small-scale chaotic flow. Manasseb^^ pointed out that the most violent 
resonant collapse occurs for the resonant Kelvin mode with the smallest radial wave number 
(called type A breakdown), leading to a fine-scale turbulence. The resonant collapse can 
display a cyclic dynamics with a transient chaotic phase followed by a relaminarization of 
the fiow^>^. During the chaotic phase, the kinetic energy remains mostly in modes with 
frequencies in the range of those of the inertial modes^S. 

It is tempting to speculate that the turbulence appears progressively from the chaotic 
set associated with the resonant collapse, but the experimental observations reports a more 
complex road to turbulence. Indeed, by further increasing the precession ratio, Lin et a!^ 
observed a non-linear regime after the resonant collapse, characterized by the disappearance 
of the free Kelvin modes in the power-spectrum. At critical precession ratio, a sustained 
intense turbulent flow fills the vessel2ii2I. Gans^S reported a discontinuity of the torque curve 
at the transition to turbulence. The torque curves exhibits also a hysteresis between the 
turbulent and non-turbulent branches. A similar bistability has been reported by Mouhali et 
in a cylinder and by Malkus^ in a ellipsoidal geometry. The presence of a discontinuity 
and the hysteresis suggests that a subcritical instability triggers the turbulence for these 
experiments with large Reynolds number {Re > 10®). The nature of the transition to 
turbulence remains unclear. 

In the framework of dynamical systems, the so-called subcriticality refers to a transition 
from a base state to a new state (here turbulent), such that the new state is created from one 
or a sequence of bifurcations^^, which does not change the local stability of the base state. It 
is analogue to a first-order phase transition, like the vapour-liquid transition of water below 
the critical point. The transition is generally abrupt and violent and the triggering depends 
on the strength of the perturbation if the flow remains linearly stable. In the past fifty 
years, great progress has been made to understand the subcritical bifurcation to turbulence 
in shear flows, like pipe fiows^^ or Plane-Couette fiows^^^— . Recently, those concepts have 
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been extended to magnetohydrodynamics . Subcritical bifurcation in shear flows 

is characterized by specihc properties^iiSi^ such as non-normal growth, complex structures 
of the flow pattern and hnite lifetimes of the turbulence^^— . 

The hnite lifetime of the turbulence is one of the most striking features of shear hows 
undergoing a subcritical bifurcation to turbulence. The perturbations leading to turbulent 
spots may relaminarize after a hnite duration. The lifetimes of turbulence are random 
and their mean value increases rapidly with the Reynolds numbers. It turns out that the 
turbulence is a metastable regime for a range of Reynolds numbers. This feature has been 
observed in hydrodynamic shear hows (pipe how^, Plane-Couette how^^ or stable Taylor- 
Couette how^) and in magnetohydrodynamics shear-hows^^*^. 

In the present study, we investigate the subcritical transition to turbulence in a huid 
how driven by a precessing cylindrical vessel. The observed transition to turbulence is a 
robust process and occurs above a well dehned threshold diherent from the threshold of the 
parametric instability of the forced Kelvin mode. Below the threshold of the subcritical 
transition, the turbulence is self-sustained, describing an hysteresis bifurcation diagram. We 
show that the turbulence is metastable in the hysteresis region and can cease suddenly. 

Section m introduces the experimental set-up and the dimensionless control parameters. 
In section IIIIl the regimes appearing before the turbulence are briehy reviewed. The quan¬ 
titative features of the turbulence and the associated subcritical bifurcation are presented 
in section IIVI Section |V] presents the experimental investigation of the hnite lifetime of the 
turbulence. 


II. THE EXPERIMENTAL DEVICE 
A. Experimental setup 

The experimental setup is illustrated on the Fig. [1] The vessel is a cylinder of radius 
R = 163 mm and height h = 326 mm. It is completely hlled with water. For the visualization 
of the how (see Figl2]), a small amount of air is introduced. The bubbles are only used to 
visualize the how pattern and all the quantitative results have been obtained with a very 
small amount of bubbles. The container rotates around its symmetry axis. The cylinder is 
driven by an asynchronous 3 kW motor via a transmission chain. The power is supplied by a 
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slip ring. To improve the parallelism of the end caps of the cylinder, they are joined together 
by eight rods (yellow rods on the right Fig{T]). The end-cap located at the opposite side of 
the motor is maintained by a counter bearing, in order to avoid any bending. This structure 
is mounted on the turntable. It can be tilted to vary the inclination between the direction 
of the rotation of the cylinder and the direction of the precession. In the present study, 
the angle is constant at 90° so that rotation axis and precession axis are orthogonal like in 
refs^Si^. The turntable is driven by a second asynchronous 2.2 kW motor. The frequency 
rate of the turntable ft can be increased up to IHz. We have not used a control loop for 
the motors, in order to avoid perturbations due to the control process, like self excitations, 
which stress the mechanics of the experiment. The rotation rates of the cylinder and the 
turntable are continuously measured by two tachometers. The fluctuations of the rotation 
rates are smaller than 0.05% and the rotation rates do not vary during the transitions from 
and to the turbulent regime. 



FIG. 1. Left: photography of the precessing cylinder. Center: transversal photography. The white 
pressure sensor is located at the bottom of the end-cap, close to the lateral-wall (white arrow). An 
array of UDV sensors is also visible. Right : sketch of the experiment. 


B. Measurement techniques 

1. Power consumption measurement 

The power consumption of the three-phase asynchronous motors is calculated by measur¬ 
ing the input current on a single phase and the voltage difference between two phases. Here, 
we report only the power consumption of the motor of the cylinder. The measurement of the 
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current is performed with a Tektronix current-clamp TCP312 and amplified by a TCPA300 
amplifier with a typical accuracy of 1%. The voltage is measured by a differential probe 
Sapphire SI-9010. The amplitudes of the current I{t), the amplitude of the voltage differ¬ 
ence U{t) and their relative phase 0(t) are extracted by a Hilbert transform. The internal 
resistance of the motor Rj has been measured by a 4 wires technique after a long run, in 
order to measure Rj close to the operating conditions (effect of the internal heating). The 
functions U{t), I{t) and 0(t) vary on time scales much larger than the rotation period of the 
motor. The instantaneous power consumption P{t) of an equilibrated motor is defined by 

Pit) = ^U{t)I{t) cos (^(j){t) + 0 - ^RiHt)^- (1) 

Because I(t) and U(t) are the peak values, the product UI and P must be multiplied by 
the factor 1/2 to obtain the RMS values. The relative phase 0(t) has been shifted by vr/b in 
order to obtain the relative phase between the current and the voltage for the same phase. 
The average power consumption is given by Pm = T~^ Pdt, with T the duration of the 
measurement. The value Pm can be decomposed into the sum of the mechanical loss, the 
magnetic loss and the power dissipated by the flow. Only the latter is of interest to us. 
We expect that the magnetic and mechanical losses do not change signihcantly for constant 
rotation rate Oc- Their values have been measured for runs at 0* = 0 (no precession), 
which corresponds to a solid body rotation without power dissipated by the flow. Our 
measurements show that the power consumption increases like for 0^ = 0. 

2. Pressure measurement 

The pressure is measured at the end-cap close to the motor (see center FigH]). Two XPM5 
miniature pressure sensors are mounted flush on the rotating end-cap. They are on the same 
radius with respect to the axis of revolution, located at r = 160 mm. Their diameter is 3.6 
mm. In the present study, we focus only on the pressure variation. In section IIIIl we present 
the frequency power spectrum. The power spectra are rescaled in order to have a power 
spectral density PSD[p]{f) tending to 1 for low frequencies. In the rest of the paper (section 
El an dCYD, we have used a smoothing hlter to low-pass hlter the signal with a typical cut-off 
frequency /c/2. The purpose of the hlter is to suppress the components with frequencies 
equal or larger than the frequency rate of the cylinder fc in order to diagnose the variation 
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of the slow components. 


C. Dimensionless number 

We recall all the parameters characterizing our set-up in the table I. 


Parameters 

Definition 

Value 

1/ 


Kinematic viscosity 

10“®m^ s“^ 

R 


Radius 

163 mm 

h 


height 

326 mm 

fc 

= rtcl2'K 

Rotation frequency 

0 - lOHz 

ft 

= Ri/27r 

Precession frequency 

0 — IHz 

0 


Precession angle 

90° 


TABLE I. Physical parameters of the experimental setup. 


The dimensional analysis allows us to dehne three dimensionless numbers from the hve 
parameters (h, R, u, ^t) (we consider 0 as a dimensionless number) 


r = A 

2R' 




e = 


Re = 




( 2 ) 


r is the aspect ratio and is equal to 1. It plays an important role in the determination of the 
eigen-frequency of the Kelvin modest. The Reynolds numbers Re quantihes the dissipative 
process and can be varied in the range [1.6 — 16.5] x 10®. The effects of the precession are 
quantihed by the precession ratio e, which varies in the range 0 — 0.16. Only the latter two 
dimensionless numbers can be modihed. In the present study, the Reynolds number is only 
changed for studying the scaling of the power dissipated at constant e (section IIV A|) and 
the parameter space (section lIV CD . In the most part of the paper, the Reynolds number is 
hxed to Re = 5.65 x 10®. 


III. THE FLOW REGIMES 

Before focusing on the transition to turbulence, we present the different regimes observed 
in the precessing experiment. We distinguish three regimes by increasing progressively the 
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precession ratio e for a constant Reynolds number Re, thanks to a direct visualization 
of the flows (see Figj2]) or by analysing the power spectrum of the pressure signal (see 
FigJS])- Visualization and pressure measurements have been performed for two different 
Reynolds numbers: Re = 14 x 10® for the visualization and Re = 5.65 x 10® for the pressure 
measurements. 

Our classihcation holds for the range Re G [1,14] x 10® and the parametric study of the 
regimes is reported in the section HV Cl (cf Fig. [7]). The three regimes are 

1. the laminar regime, 

2. the non-linear regime, 

3. the turbulent regime. 



(a) (b) (c) (d) 


FIG. 2. The different regimes of the precessing flow are observable, via a small amount of bubbles, 
in the precessing reference frame (turntable frame) for R = lOHz and e = 0.74x 10“^ (a), 5.9 x 10“^ 
(b), 7.4 X 10“^ (c), and 10 x 10“^ (d). The amount of air is kept constant. From left to right : 
standing forced Kelvin mode with a S-shape (a), non-linear regime with almost not distinguishable 
Kelvin modes and the bubbles localized in the upper left region (b), turbulent regime with an 
important spreading of the bubbles (c), strong vortex superimposed with turbulent fluctuations 
(d). 


The laminar flow corresponds to a superposition of Kelvin modes with an azimuthal 
wave number m = 1. Following the theory of precessing flow^, the hrst resonance occurs 
at / = 0.996/c for an aspect ratio h = 1. Our conhguration corresponds to a resonant case 
with one dominant Kelvin mode. In the Fig. Ha), we clearly identify the S-shape of the 
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minimal pressure region, associated with the axial wave number kz = n ■ L~^. It is standing 
in the turntable frame. The power spectrum of the pressure signal exhibits an intense peak 
at the frequency /c and at the corresponding harmonics (see FigJ31(a)). The peaks at 2/c 
and 3/c likely arise from non-linear self-interaction of the forced Kelvin mode^i. For a forced 
Kelvin mode defined by the frequency /c , an azimuthal wave number m = 1 and an axial 
wave number ko , the non-linear interaction will drive the modes (2/c, m = 2, fc = 0) and 
(0, 0, fc = 2/co) (the other combinations are forbidden). Moreover, the inertial pressure term, 
given by |r2 x r|^, contributes to the component / = 2fJ^. The centrifugal pressure and the 
hydrostatic pressure contribute also to the pressure components at fc- 



FIG. 3. Power spectrum of the pressure signal for Re = 5.65 x 10® and e = 0.68 x 10“^(b) 
(laminar regime), e = 5.47 x 10“^(b) (non-linear regime with triadic resonance), e = 6.84 x 10“^(c) 
(non-linear regime without triadic resonance) and e = 9.57 x 10“^ (d)(turbulent regime). 


We dehne as non-linear regime the regime between the hrst instability of the Kelvin mode 
m = 1 and the transition to turbulence. At a given precessing ratio, the forced Kelvin mode 
becomes unstable. Indeed, the peaks observed on the power spectrum (see Figj3](b)) suggest 
the existence of two free Kelvin modes (/i = 0.223/c and /2 = 0.75/c) in near-resonance 
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with the forced Kelvin mode (/i + /2 = 0.97/c). The thresholds of appearance of these 
triads will be discussed later (Fig. [7]). At the threshold, the peaks in the power-spectrum 
correspond to frequencies close to /i = 0.35/c and /2 = 0.65/c, which may be assigned to 
m = 5 and m = 6 modes. By further increase of the precession ratio, these frequencies vary. 
We believe that this effect is due to the saturation process which is achieved by a detuning 
in frequency of the modes (see for instancei^). These results will be reported in a future 
paper. 

Further increasing the precession ratio, we have identihed a sequence of bifurcations, and 
we are currently designing a velocimetry set-up to diagnose in details these bifurcations. 
These results will be reported in a future paper. For larger precession ratio, the two previous 
frequency peaks disappear (see Figj3](c)) and the power-spectrum becomes flatter. The 
forced Kelvin mode is no more distinguishable by direct visualization (see Figj2](b)). The 
bubbles are localized close to one end-cap due to the breaking of the centro-symmetry of 
the laminar flow. 

These observations suggest that the dynamics is already complex and likely chaotic. 
Our measurements are consistent with the observations of Lin et alM, who observed the 
disappearance of triadic resonance of inertial modes in the strongly non-linear regime. 

Further increasing the precession ratio, the flow becomes turbulent. The properties of 
this transition will be detailed in the following sections. Visually, the turbulent flow is 
characterized by the spreading of small bubbles in the full vessel (see Fig|2](c)). Occasionally, 
we observe a vortex hlament aligned with the direction of precession (see FigJ21(d)). The 
pressure signal gives few information about the turbulent flow (see Fig. [3](d)). The power- 
spectrum becomes almost flat on the range of frequencies / < 4/c. 


IV. TURBULENCE AND SUBCRITICAL BIFURCATION 

In this section, we focus on the transition to turbulence. First, we conhrm the existence 
of a discontinuity and a hysteresis of the power consumption at the transitiort^^. Then, we 
study the dynamics and the robustness of this transition. We conclude this section by the 
systematic determination of the threshold in function of the Reynolds number. 
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A. Power consumption: hysteresis and scaling 


We carried out measurements of power consumption at constant Reynolds number Re = 
5.65 X 10^ (/c = 4Hz). We have performed a first ascending series of measurements from 
e = 0 to 16.4 X 10“^ and a second descending series from 10 x 10“^ to 6.56 x 10“^. Due to the 
presence of fluctuations at low frequencies (see Figj3]), each measurement lasts 40 minutes 
(9600 periods) in order to converge the average power consumption Pm{^)- The results are 
reported on Fig. 01 
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FIG. 4. Mechanical power consumption of the motor of the rotating cylinder for Dc = (Svr) rad/s 
and e = Rt/Dc € [0,0.165]. Two distinct branches are clearly observable: the lower branch (L) 
with the laminar regime (circle) and the non-linear regimes (square), and the turbulent branch 
(r) (triangle). Dashed line: linear fit with Pm = PmO + Ce. Insert: zoom on the hysteresis. The 
system jumps from the L branch to the T branch at elt = (7-3 ± 0.1) x 10“^ and goes back from 
T to L at etl = (6.7 ± 0.1) x lO’T 


The average power consumption Pm{^) displays two distinct branches: the lower branch 
L (squares and circles) with e G [0, and the turbulent branch T (triangles) with 
e G [ctl, 16 X 10“^]. The L branch corresponds to the laminar regime (circle) and to 
the non-linear regime (square). Both branches are disconnected and the system undergoes a 
subcritical bifurcation from L to T at and goes back from T to L at e^L. The thresholds 
are e^T = 7.3 ± 0.1 x 10“^ and e^L = 6.7 ± 0.1 x 10“^ and the hysteresis is clearly visible 
with clt — ^TL — 6 X 10 “^ (insert, Fig0]). The determination of eir is very robust (section 
IIVB|) whereas a precise measurement of e^L is very difficult [section IV A|) . 
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The behaviour of the power consumption of the L and T branches are quite different. On 
the L branch, the power consumption does not vary signihcantly for e < 0.05 and increases 
rapidly between e ~ 0.05 and e^T- In this range, the resonant triads disappear (section 
HID. The power consumption on the T branch exhibits a linear growth (dashed grey line on 
FigSD with 


P^{e)=Po + Ce (3) 

with Pq = SOW and C = 124hT. A similar behaviour was observed in the torque measure¬ 
ment of Gans^S for the hydrodynamic case. 

We performed another series of measurements with constant e = 8.75 x 10“^ in order to 
estimate the Reynolds dependence of the power consumption by varying The flow re¬ 
mains turbulent for the associated Reynolds numbers, which vary in the range [3.5,13] x 10®. 
To compare the measurements, we subtract the power consumption at e = 0 to remove 
the mechanical loss (increasing like of the set-up. The power injected in the flow 
AP = Pm(e) — Pm(0) is reported on the left hand side of Fig. |3 It increases like G® 
for constant e. The dimensionless power AP/Pg, with P* = is roughly con¬ 

stant for Re G [3.5,13] x 10® (right FigJSD, suggesting that the dimensionless injected power 
becomes independent of the Reynolds number. 
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FIG. 5. Left: power injection, AP = Pm(e) — Pm(0) for e = 8.7 x 10 ^ and Re E [3.5,12] x 10®. 
The dashed grey curve follows the power-law AP ~ Rj?. Right: compensated power AP/Pg with 
Pg = p2TTR^Lnl. 
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From equation [3] and the behaviour shown in Fig. [5l it can be concluded that the power 
dissipated into the turbulent flow scales like 


AP 


pR^mi ^ ^ 

This scaling differs from the linear viscous theory with AP/Pg ~ e^Re and the theoret¬ 
ical constant bounds, AP/Pg ~ e°, expected for large Reynolds numbers and an infinite 
cylinder—. Additional information about the turbulent flow is required to explain the 
e-dependence of the above scaling law by dimensional analysis. 


B. Transition to the turbulence 

We have performed a series of measurements to investigate the transition to turbulence 
at the threshold at a given Reynolds numbers Re = 5.65 x 10^. The precession ratio is 
increased step by step with a plateau time of 60 seconds (240 f~^) and a precession ratio step 
Ae = 1.42 X 10 We performed 16 runs with the same procedure, using the pressure signal 
for the detection of the transition (see left hand side Figj6]). The transition always occurs 
when the precession ratio reaches the value e = 7.25 x 10“^. The black curve represents 
the precession ratio (right ordinate) and the blue curve and the red curve correspond to the 
pressure signal from the same probe but for two different runs. The pressure can slightly 
drift due to the thermal expansion of the water but we did not observe any correlation 
between the drift of the pressure and the onset of the instability, so we believe that this 
effect does not play any significant role in the dynamics of the flow. We observe at td = 294s 
(resp. td = 309s) for the blue (resp. red) curve a drop of the pressure after the growth of 
the precession ratio to e = 7.25 x 10“^ at tp = 280s. It corresponds to the threshold e^T 
for the transition to the turbulence. We have measured the delay AT = td — tp between 
the change of e and the beginning of the instability. The average value of the delay AT is 
15 seconds (60 fp^) with a standard deviation of 9 seconds (0.6 x AT). The delay AT is 
comparable to the precession period = 3.3s and could correspond to the time needed 
to re-establish the base flow, which then becomes unstable. In any case, the delay AT is 
much smaller than the typical transient life time of the turbulence in the hysteresis region 
(see section lYD- 

The discontinuity of the pressure confirms the presence of a subcritical transition. The 
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pressure jump corresponds roughly to Ap ~ 500 Pa and is detected simultaneously on 
the diametrically opposed probes. The transition seems to be global and affects the low 
frequency components of the pressure field. The pressure drop is likely due to a slowdown of 
the rotation of the bulk flow. An estimation of this slowdown gives A/ = (Ap/ {2'kR) ~ 
0.7 Hz, corresponding to 17.5% of the imposed frequency rate fc- The slowdown of the bulk 
flow in the turbulent regime has also been reported for processing cylinders^! and spheres^S. 



FIG. 6. Left: evolution of the low-pass filtered pressure signal p (blue and red) for two different 
realizations, and precession ratio e (black) as a function of time t. Right: rescaled pressure tt (equa¬ 
tion [5]) during the transition. The different realizations are in grey, and the averaged realization is 
in black. 


In order to compare the different realizations, we define a rescaled pressure tt such that 


7l{t) = 


p{t) - (pt) 


(5) 


(Pl) - (pt) 

with {pl) (resp. (pt)) the average pressure before (resp. after) the transition to turbulence. 
The average pressure vr decreases from 1 to 0. We have used a coherent average to super¬ 
impose the curves. It consists in finding the maximum of the correlation function of two 
different realizations, in order to remove the delay between them. The instant t = Os cor¬ 
responds to the beginning of the transition. The different realizations (right hand side Fig. 
El grey curves) are concentrated around the average curve (black curve) and do not display 
any exponential growth or exponential relaxation to the turbulent branch. We estimate the 
duration of the transition Tj = 4 ± 0.5s (16 /“^). 
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C. Effects of the Reynolds numbers 


Now, we investigate the dependence of the threshold as a function of the Reynolds num¬ 
bers Re . The covered parameter space {Re, e) is reported on the left hand side of Fig. [3, 
with eiT (circle) and ctl (square). The diamonds represent the frontier between the laminar 
flow and the non-linear regime. It corresponds to the appearance of free Kelvin modes in the 
power spectrum fsection HTB . The region between e^T and ctl corresponds to the bistable 
region between the L and the T branch. The threshold e^T depends weakly on the Reynolds 
number, following an empirical power law with eir Re (Fig. [3 right). The width of 
the hysteresis region given by the difference e^T — ^tl, is roughly independent of Re. 
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FIG. 7. Left : parameters space Re versus the measured thresholds of the first instability of the 
forced Kelvin mode (diamonds), the threshold elt of the transition from L branch to T branch 
(circles) and etl from T to L (squares). Right : zoom on the behaviour of elt and etl- Dashed 
curves show a fit to a power law with e ~ . 


V. FINITE LIFE TIME OF TURBULENCE 
A. Metastability 

The turbulent flow exhibits an interesting feature on the range of the hysteresis e G 
R can suddenly cease and relax to the L branch, with a random lifetime in 
the turbulent branch. Figure [8] illustrates this phenomenon. The time series of the power 
consumption for three different runs at Re = 5.65 x 10^ are displayed on the left hand 
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side. For t < 0 s, the flow is turbulent with e = 9.3 x 10“^ (> eLr) and at t = 0 s, the 
rotation rate of the turntable is reduced to e = 6.84 x 10“^ with ctl < e < eLT- The drop 
of power consumption corresponds to the transition from the T branch to the L branch. 
It occurs rapidly after the slowdown for the red curve, or 400s (1.6 x 10^/“^) for the blue 
curve. The third realization (black curve) is still turbulent at t = 700 s (2.8 x 10^/“^). 
Unlike the resonant collapses'll, the relaxation to the L branch is definitive. We did not 
observe spontaneous transition from the non-linear regime to the turbulent regime below 
the threshold elt- 



FIG. 8. Three different realizations (red, blue, black curves) of the relaxation from the T to L 
branches for Re = 5.65 x 10^ and e = 8.5 x 10“^ (> clt) for t < Os and e = 6.84 x 10“^ (< clt) 
for t > Os. Left: AP(t) power consumption. Right: rescaled pressure 7r(t). 


The pressure measurements display also an abrupt transition when the turbulence col¬ 
lapses. The corresponding data are reported on the right side of Fig. [HI The rescaled 
pressure (defined in section IIVBD suddenly increases at the transition. The dynamics of this 
relaxation is discussed in section IV Cl The turbulence seems to be metastable for a given e in 
the hysteresis region [e^L, glt]- The lifetimes of the turbulence are random and a statistical 
study of the lifetimes is performed in the next section. 

As mentioned in the section lYl the precise determination of e^L becomes intractable due 
to the apparent metastability of the hysteresis part of the turbulent branch. We carried out 
measurements for e = 6 x 10“^,i.e. below ctl, with (e — eTL)/^TL — 10%. The transition 
is almost immediate. The delay lasts 10 s (40 /t^) with a small standard deviation of 1.5 
s. This time scale is comparable to the delay AT before the beginning of the transition 
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to turbulence fsection lIVBD . It turns out that the turbulence is no more metastable for e 
smaller than ctl- 


B. Lifetime statistics 
1. Lifetime distribution 


To quantify this stochastic process, we have investigated the statistical properties of the 
duration of the turbulence in the metastable regime. The experimental procedure consists 
in starting a run at a precession ratio e = 9.3 x 10“^ (> slt = 7.25 x 10“^) in order to 
obtain a strong turbulent flow. After 200s (8 x 10^/“^), the rotation rate of the turntable 
is reduced to the e wanted. The lifetime r is dehned as the duration between the end of the 
slowdown and the pressure rise, corresponding to the relaxation to the L branch (see hg[8]). 
The incertitude concerning the end of the slowdown is of order AT (cf. llVj) . 

We have performed a large number of measurements (N = 125) for e = 6.7 x 10“^ in 
order to calculate the probability distribution function p(r) (left Fig. [9]). The distribution 
displays an exponential tail with p(r) ~ exp(— t/tq) and tq = 288±30 s or 1.1 x 10^ rotation 
rate periods (left Fig. [9]), which is characteristic of a memoryless process. The value Tq is 
close to the average lifetime (r) = 267 s, directly calculated from the set of durations r. 
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FIG. 9. Left: distribution function of the lifetime of the turbulence for e = 6.7 x 10“^ and 
Re = 5.65 X 10^ with P{t) = Cexp(— t/tq) with tq = 288 ± 30s. Each point corresponds to an 
interval of 150 s. Right: ratio of the standard deviation cJt- over the average (r). 
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2. Average lifetime 


We have performed a series of measurements to calculate the average lifetime in the 
turbulent branch for e between 6.67 x 10“^ and 6.84 x 10“^ (Fig- [lO] )• The average is 
calculated for at least 30 — 40 realizations for each precession ratio. The average lifetime (r) 
varies from 170 to 2110 seconds. It exhibits a steep increase when e comes closer to e^T- 
By assuming that (r) becomes inhnite when e is equal to elti it is tempting to model the 
behaviour of the average lifetime (r) by a critical behaviour with ~ |e — eLT\~^■ Here, 
the critical point is ”a priori” dehned as equal to elt but it cannot be ruled out that the 
critical point differs from this value. The exponent (3 estimated from such a £t turns out to 
be larger than 10, suggesting that the divergence of (r) is stronger than a power law. The 
behaviour of (r) is better approximated by an exponential growth with 

= exp(a|e - + fe) (6) 

and (a, b, (3) = (1.08, —9.13,1/2). The £t is represented by a red curve on the left hand side 
of Fig. [TUI Equation (|U]) is similar to the Arrhenius equation and may refer to a behaviour 
called ’’supertransient”— in dynamical systems. 
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FIG. 10. Left: Average lifetime of the transient turbulence for Re = 5.65 x 10® for different e < slt 
{ eLT = 7.3 X 10“^). The red curve represents the functional form (r)t/i = exp(a|e — + b) 

and the blue curve corresponds to the exponential function {T)th = exp(ae + b). Right: Average 
lifetime (r) in function of |e — 
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We have represented the escape rate (r) as a function of |e — eLT\~^'^ in a logarithm 
scale for the ordinate (right Fig JTOj) . The circles are aligned along the red line, justifying the 
functional form chosen. The coefficients a and b have been estimated from the linear relation 
log((r)) = a|e — + b with a least square method, which gives a smaller relative error 

(cf. equation [7]) . 

Actually, the value (3 can be varied in the range [0.1, 0.7] without changing signihcantly 
the quality of the fit. We define the relative error function re as 


{r)l ) 


(7) 


with N = 7 the number of experimental points. The relative error Ve is between 9.45% and 
9.55% for (3 G [0.1, 0.7] and all the functions {T)th collapse almost on the same curve. The 
main contribution to the relative error r comes from the mean durations (r) at e larger than 
6.8 X 10“^. In order to calculate more precisely the exponent /3, measurements at larger e, 
i.e. with much larger (r), would be needed. The quality of the £t depends also on the value 
of the critical point, here dehned ”a priori” as e^T- 

We point out that the behaviour of (r) can be also approximated by an exponential 
function with {T)th ~ exp(ae) with a relative error = 7% (Right FigHOl blue dashed line). 
The £t is even better than with the equation El but it implies that the turbulence remains 
always metastable. 

The difference between the relative error of the exponential law and the equation El comes 
from the contribution of the two last points, where the incertitude on the measured (r) is 
large (around 10%). It implies that our limited set of measurements does not allow us to 
decide if there is a threshold (cf equation E]) or if the turbulence remains metastable for larger 
e (exponential law). An estimation of (r) for eiT with the exponential law gives (r) = 10® s 
or 11 days. It implies that the experimental investigation of (r) near slt is impossible with 
our set-up. 

With the same set of measurements, we have calculated the standard deviation of (r) 
for different e. The standard deviation normalized by (r) (right Figj9]) is close to one, as 
expected for an exponential distribution. 
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C. Relaxation process 


The relaxation from the turbulent to the lower branch is a fast process, compared to 
the transient lifetime of turbulence. We noticed by direct visualization that the turbulent 
flow still fills the entire vessel (like in the right hand side of Figj2]) before the collapse of 
the turbulence. We have reported different realizations (in grey) of the relaxation of the 
pressure tt on the left of Fig. [Tl]for Re = 5.65 x 10^ and e = 6.7 x 10“^. The procedure is 
identical to the one detailed in section IIVBI 
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FIG. 11. Left : superposition of different realizations (in grey) of the relaxation process for Re = 
5.65x10® and e = 6.7 X 10 ^ of the low-pass filtered pressure. The black curve is the average, which 
follows an exponential function (insert). Right: the relaxation time 7 “^ versus the precession ratio. 


All the realizations are concentrated around the average (left hand side of FigJTTl black 
curve). Unlike the stochastic process of the lifetimes, the relaxation seems to be deterministic 
with a small variability in the realizations. The transition lasts roughly 20 seconds. The end 
of the process {t > 5s) follows an exponential relaxation with tt — 1 ~ exp(— yf) and 7 > 0 
(insert in left hand side Fig JTT]) . The relaxation times 7 “^ are reported for different e on 
the right of Fig. [TT] The parameter 7 “^ is almost constant and stays in the range [6.5, 7.5] 
seconds. For e = 6 x 10“^ (< eTL)i fhe relaxation remains also exponential with a similar 
time scale 7 “^ = 7.69s, even if the turbulent state is no more metastable in this range of e 
(cf section rvWI) . 
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VI. CONCLUSION 


A. Summary 

In the present study of a precession driven flow, we have conhrmed the existence of 
a subcritical bifurcation to turbulence associated with a discontinuity and a hysteresis. 
The observed transition to turbulence is a robust process and occurs above a well defined 
threshold eLT, different from the one of the parametric instability of the forced Kelvin mode. 
Both the pressure and power consumption measurements exhibit an abrupt variation at the 
transition from and to turbulence. The behaviour of the power consumption changes at 
the threshold, with a linear dependence of the power on the precession ratio. The pressure 
measurements suggest that the transition is associated with a slowdown of the solid body 
rotation. The threshold e^T of the transition depends weakly on the Reynolds numbers. 

In the region of the hysteresis e G [ctliClt], the turbulence can suddenly collapse and 
lead to a definitive relaminarization of the flow. For a given precession ratio and Reynolds 
number, the lifetimes r of the turbulence are random and are exponentially distributed . 
The average lifetimes (r) of the turbulence increase rapidly when e approaches the threshold 
eLT of the sustained turbulence. The growth of the lifetimes r is well modelled either by 
the functional form (r) ~ exp(a|e — eLT\~^) with a and (3 positive or an exponential law. A 
detailed study of the relaxation process shows that the slow component of the pressure follows 
an exponential function with a relaxation rate weakly dependent of e, even for e < e^L. 

As a conclusion, a schematic bifurcation diagram is presented on the Fig. [121 Three 
regimes are observed : the laminar and the non-linear regimes on the L branch and the 
turbulent regime on the T branch (cf section |TTT|). The L branch becomes unstable at eir 
and the system undergoes a subcritical bifurcation to turbulence. The flow comes back to the 
L branch for e smaller than ctl (cf section [IV]). Between ctl and eLT, the turbulent regime is 
metastable and exhibits transient lifetimes (cf section jV]). The sneaking arrow corresponds 
to the relaxation process to the lower branch. The grey line represents the frontier between 
the basins of attraction of each branch, even if the turbulence is not strictly speaking an 
attractor. The collisions between the grey and blue lines (black dots) are only symbolic and 
do not refer to a saddle-node bifurcation. 
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Turbulent branch 



FIG. 12. Schematic bifurcation diagram of the precessing cylinder for e < 16 x 10 ^ and Re E 
[ 10 ^ 10 ®]. 


B. Discussion 


The observation of a well-defined threshold for a snstained turbulence may lead to 
different conclusions. It could indicate the presence of a linear instability past the onset 
eiT- The instability could also be triggered by non-normal and non-linear mechanisms. 
Moreover, if we consider the flow as already chaotic before the instability, the transition 
could be explained by a crisis^ of this chaotic set. This bifurcation refers to the collision 
of the basin of attraction of two chaotic attractors, leading to the emergence of only one 
attractor. However, a non-linear instability is not consistent with the measured robustness of 
the threshold, whereas the crisis would be associated with random waiting times before the 
transition. The linear instability is thus the best candidate for the transition to turbulence. 
As the instability of the laminar flow occurs for smaller precession ratios than e^T, the 
mechanism of this new instability could differ from the parametric instability of the Kelvin 
modes^. 
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